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Abstract 



I review some of the history of the flux-tube model, concentrating on work 
done by the Los Alamos-Tel Aviv collaboration on particle creation and back- 
reaction in uniform fields and in the central rapidity region. I discuss the 
^ ' incorporation of more realistic geometry and structure of the flux tube via 

application of the dual superconductor model of confinement. 



^ , During the decade of my collaboration with Judah, his main area of interest 



was the quark-gluon plasma and the efforts to create it in relativistic heavy- 
ion collisions. He was especially drawn to the problem of particle creation 

■ in the earliest stages of a collision, and the process by which these particles 
, might reach thermodynamic equilibrium. For a model of particle creation, we 

0^ ' focused on the flux-tube model. It presented us with the practical question of 

- connecting fundamental quantum field theory with phenomenology, and thus 

with the theoretical challenge of connecting quantum field theory with kinetic 
Oh, theory and hydrodynamics.^ 

r-| ! 1 Kinetic theory of particle creation 

. ^ ^ After years of effort and anticipation, the Relativistic Heavy Ion Collider (RHIC) 

■ at Brookhaven will soon begin producing collisions of gold nuclei in colliding 
H \ beams, at an energy of 100 GeV per nucleon in each beam. Theoretical con- 

siderations, as well as data from the successful heavy-ion program at CERN's 
SPS, lead us to expect that in a central Au-Au collision the incoming baryons 

^Invited talk at Nuclear Matter, Hot and Cold: a Symposium in Memory of Judah M. 
Eisenberg, Tel Aviv, Israel, April 14-16, 1999. Work supported by the Israel Science Foun- 
dation under Grant No. 255/96-1 and by the Basic Research Fund of Tel Aviv University. 

^We enjoyed a long-term collaboration in this area with Fred Cooper and Emil Mottola of 
Los Alamos and with Yuval Kluger, who wrote his Ph. D. thesis on the subject at Tel Aviv 
and then moved to Los Alamos. 
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will pass through each other with minimal deflection. The net outgoing baryon 
number will thus be concentrated within a couple of units of rapidity of the 
beam rapidities, leaving a large central region to be filled by a baryon-free 
plasma. Energy deposition in the central region will take place mainly via soft 
mechanisms, although minijet production will play a significant role. 

A QCD-inspired model for soft production in the central region begins [Q 
with the exchange of soft gluons between the two nuclei as they pass through 
each other. This leaves them charged with color, sources of color electric flux. 
The flux fills a tube connecting the receding nuclei, which thus form a "color 
capacitor" if the field is coherent across the tube, or a "color rope" if it is not. 
The electric field strength in either case is proportional to the square root of 
the number of gluon exchanges per unit area, and thus E ^ A^^^, as a result 
of a random walk in color space. 

An electric field will pop pairs of charged particles (in this case, colored 
particles) out of the vacuum, especially if they are light on the scale of the 
available field energy density. Casher, Neuberger, and Nussinov |^ derived via 
WKB the rate for particle creation per unit volume as a function of transverse 
momentum. 



= e£,log 



dtdVcPp 



1 + cxp ^ -i^ 

eh 



(1) 



(Integrating (|^) over p± gives the famed Schwinger formula [^.) 

If the flux tube is narrow to begin with, as in the case of e"'"e~ annihila- 
tion into quarks or of pp scattering, the quarks produced via (|^) will break 
the tube and begin an inside-outside cascade. For nucleus-nucleus collisions, 
however, the geometry more closely resembles that of the idealized problem of 
a homogeneous electric field created by infinite, parallel capacitor plates. Then 
the created particles form a uniform current that screens the field gradually 
according to (Abelian) Maxwell's equations. 

This is called back-reaction. A straightforward approach to the particle dy- 
namics is given by kinetic theory, where the pair creation rate (^ appears as a 
source term for the Vlasov equation M p| , 



df df 

-£ + eE-^ ^ eE\og 

ot apz 



1 -I- exp 



eE 



&{Pz) ■ (3) 



(The particle density / feeds back into the Maxwell equation (g) through the 
current.) This gives a beautiful description of matter creation and subsequent 
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plasma oscillations. The drawback of this approach, of course, is that the use 
of the classical Vlasov formalism lacks fundamental justification. 



2 Enter quantum field theory 

Quantum field theory allows a first-principles approach to the problem ^. 
One begins by taking the number of flavors Nf to infinity in order to justify a 
classical, mean-field approximation for E. If we begin with a scalar field in one 
dimension, the ingredients of the theory are then the Klein-Gordon equation 
for the matter field, 



and the semiclassical Maxwell equation, 

d.F^" = (0|r|0) . 

A Fourier expansion of 

dk r 



27r 



fkit)a, + r,{t)bl. 



ikx 



inserted into 



where 



gives evolution equations for the Fourier amplitudes fk. 



dt^ 



ojk(tf = [k - eA{t)f + . 
The canonical commutation relations for $ imply that fk can be written as 



fk{t) 



^-i J* nk{t'}dt' 



Inserting this into (0) gives equations for the effective frequencies 



(4) 
(5) 

(6) 

(7) 
(8) 

(9) 



3 

4 [nk 



(10) 



We can write the current j in terms of the amplitudes fk, and (|^) becomes the 
subtracted, renormalized equation 



A 



dk 
2^ 



{k - eA) 



(11) 
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Equations ( p^ ) and (|T^) can be integrated numerically. 

We applied this formalism to the creation of bosons ||^ and fermions Q in 
one and three Q dimensions. For illustration, I show in fig. ^ the time evolu- 
tion (in scaled variables) of a system of bosons in three dimensions. Particles 
created at early times are accelerated into plasma oscillations, with continued 
particle creation whenever the field is non-zero. A comparison to the results of 
kinetic theory [a slight modification of (||): dashed curves] shows that the lat- 
ter approximates the field-theoretic solution well. Space limitations prevent me 
from showing the beautiful results of comparing the phase space distributions 
of produced particles in the two approaches. 
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Figure 1: Plasma oscillations of bosons in three dimensions 



4 



3 Extensions 



The idealized problem of infinite, static capacitor plates was of course only a 
first step in development of the flux tube model. With an eye towards nucleus- 
nucleus collisions, we extended the analysis to the problem of capacitor 
plates separating at the speed of light, z — ±ct. As is natural in the context 
of the central rapidity region we assumed that the initial conditions are 
invariant under longitudinal boosts so that the time development depnds only 
on the invariant, comoving time t = y/t"^ — Z'^. The phase-space distribution 
of the produced particles now gives the rapidity distribution dN /dy, and anal- 
ysis of the energy-momentum tensor gives the energy density and equivalent 
temperature as functions of r (though there is as yet no mechanism for equi- 
libration). As it turns out, the results are again well approximated by kinetic 
theory. 

The geometry to this point is still that of an idealized system of infinite 
transverse extent. Upon solving the problem in a cylinder of finite radius (a 
computationally formidable problem due to nonlinear mode mixing), Judah 
found 1 12 that the agreement between field theory and kinetic theory is at best 
qualitative. 

As mentioned, the semiclassical approximation presented above finds its jus- 
tification in the limit of large Nf. The next order in 1/A''/ introduces particle- 
particle collisions and thus a mechanism for relaxation to thermal equilibrium 

. Unfortunately, the resulting integro-differential equations are (so far) com- 
putationally unmanageable. Kluger, Mottola, and Eisenberg Q later eschewed 
this approach and showed, in what was to be Judah's last paper, that equili- 
bration can arise even in the lowest-order mean field theory shown above if one 
allows for dephasing of the highest-frequency quantum oscillations. The 
connection between reversible quantum mechanics and irreversible kinetic the- 
ory has long been an area of mystery and of hard work. In this last paper, 
a Vlasov equation non-local in time is actually derived for the pair-creation 
problem]^ 



4 The flux tube of the dual superconductor 

Judah, as noted, applied the field-theoretic formalism to particle production in 
a cylindrical fiux tube of finite, fixed radius. The QCD flux tube, however, is a 
dynamical object, governed by field equations, 't Hooft and Mandelstam noted 

■^The paper is a mathematical tour de force of Airy functions, uniform asymptotic expan- 
sions, etc. Anyone who knew Judah's great love for mathematical physics of the Morse & 
Feshbach variety will recognize Judah walking in its pages. 
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long ago that flux tubes arise naturally in superconductors via the Meissner 
efl'ect; just as magnetic monopoles would be confined by a magnetic flux tube 
in a superconductor, quarks with their electric color charge would be confined 
by an electric fiux tube if the QCD vacuum has the structure of a dual super- 
conductor. Melissa Lampert and I have taken the first step of studying the 
dynamics of classical charges moving in an electric fiux tube and the reaction 
of the flux tube in this model [ p^ . 

To specify the dual superconductor model, we begin with Maxwell's equa- 
tions coupled to both magnetic and electric currents, 

d^F^" = (12) 
d^F^" = J-;. (13) 

Eq. ( p^ ) is no longer just a Bianchi identity; thus a vector potential can be 
introduced only if a new term is added to take care of the magnetic current, 

pf.^ = dt^A" - ^''A^' + e^'^'^'^Gxa , (14) 

G'^- = -n'^in ■ dr'fg . (15) 

This vector potential can be coupled to electric charges as usual; in order to 
introduce magnetic charges, one introduces a dual potential ||l^ via 

pf.^ ^ gt^^u _ ^ ^^^'^^<yMx„ , (16) 

Mf"" = -nf^in ■ d)-^f^ . (17) 

Now we can write a model for the monopoles, for which the simplest is an 
Abehan Higgs theory ||, |l9|, 

D^^^'^V + AdVP - t;')V' = , (18) 

where 

Df^ ^d^- ^gB^ . (19) 

This theory should produce the desired magnetic condensate to confine electric 
charge. The magnetic current appearing in ( p^ ) is 

^2glmi;*D''^ij . (20) 

For the electric charges, we take simple two-fiuid MHD (see for details). 

The parameters of the theory may be adjusted to put the superconductor 
into the Type I or the Type II regime; the density of the charged fluid will make 
the plasma frequency uip larger or smaller than the mass my of the gauge boson 
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Figure 2: Plasma oscillations in the Type I dual superconductor with LUp < my: 
the electric field on the axis of the flux tube. 



in this Higgs theory. I show in fig. ^ plasma oscillations in the flux tube for a 
Type I superconductor with ujp < my- The oscillations are made nonlinear by 
the reaction of the flux tube, which tries to close up when the field is weak. 

This is a first step towards the inclusion of the dynamics of the fiux tube in 
a model of particle creation. The latter might be included, as discussed above, 
as a kinetic theory based on a Vlasov equation or via a full field-theoretic 
treatment. As argued in however, the simplest Higgs theory (|l8|) doesn't 
contain enough dynamical confinement physics to produce good phenomenology. 



References 

[1] T. S. Biro, H. B. Nielsen, and J. KnoU, Nucl. Phys. B245, 449 (1984). 

[2] A. Casher, H. Neuberger, and S. Nussinov, Phys. Rev. D 20, 179 (1979); 
ibid. 21, 1966 (1980). 

[3] F. Sauter, Z. Phys. 69, 742 (1931); W. Heisenberg and H. Euler, Z. Phys. 
98, 714 (1936); J. Schwinger, Phys. Rev. 82, 664 (1951). 

[4] K. Kajantie and T. Matsui, Phys. Lett. 164B, 373 (1985). 



7 



[5] A. Bialas and W. Czyz, Phys. Rev. D 30, 2371 (1984); ibid. 31, 198 (1985); 
Z. Phys. C28, 255 (1985); Nucl. Phys. B267, 242 (1985); Acta Phys. Pol. 
B 17, 635 (1986). A. Bialas, W. Czyz, A. Dyrek, and W. Florkowski, Nucl. 
Phys. B296, 611 (1988). 

F. Cooper and E. Mottola, Phys. Rev. D 40, 456 (1989). 

Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. Lett. 67, 2427 (1991). 

Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. D 45, 4659 (1992). 

Y. Kluger, J. M. Eisenberg and B. Svetitsky, Int. J. Mod. Phys. E2, 333 
(1993). 

F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola and B. Svetitsky, Phys. 
Rev. D 48, 190 (1993). 

F. Cooper, G. Frye and E. Schonbcrg, Phys. Rev. D 11, 192 (1975); 
J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 

J. M. Eisenberg, Phys. Rev. D 51, 1938 (1995). 

F. Cooper et al, Phys. Rev. D 50, 2848 (1994). 

Y. Kluger, E. Mottola and J. M. Eisenberg, Phys. Rev. D 58, 125015 
(1998). 

S. Habib, Y. Kluger, E. Mottola and J. P. Paz, Phys. Rev. Lett. 76, 4660 
(1996). 

M. A. Lampert and B. Svetitsky, ^ep-ph/9905455| (May, 1999). 

D. Zwanziger, Phys. Rev D 3, 880 (1971); M. Blagojevic and P. Senjanovic, 
Phys. Rep. 157, 233 (1988). 

T. Suzuki, Prog. Theor. Phys. 80, 929 (1988); S. Maedan and T. Suzuki, 
ihid. 81, 229 (1989); T. Suzuki, ihid. 81, 752 (1989). 

H. Suganuma, S. Sasaki, and H. Toki, Nucl. Phys. B435, 207 (1995). 



8 



